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The role of neutrality constraints for the phase structure of QCD at non-vanishing chemical 
potentials is studied within a self-consistent truncation scheme for the Dyson-Schwinger equation 
of the quark propagator in Landau gauge. We find the (approximate) color-fiavor-locked phase to 
be energetically preferred at all potentially relevant densities and for physical values of the quark 
QQ ' masses. We furthermore observe the impossibility to define this phase by residual global symmetries 

I and discuss the role of chemical potentials. 
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The phase diagram of quantum chromodynamics (QCD) is expected to contain color-superconducting regimes at 
sufficiently high densities and low temperatures [U, [3, H, 0, [1, 0, Q • Considering three (degenerate) quark flavors, the 
fSJ I color-flavor locked (CFL) phase Q is the ground-state for vanishing temperatures and asymptotically large densities [§, 
[Toj . At densities of potential relevance for compact stellar objects it is, however, not obvious how the splitting of Fermi 
surfaces due to finite quark masses and neutrality constraints is influencing the ground-state [ll|, [l^, ■ In Ref. [l^ 
^ ^. we have elaborated the role of the renormalized strange-quark mass when addressing this question. Including medium 
I ' modifications of the effective quark interaction decreases the dynamically generated contribution to the strange-quark 
mass. Therefore, for the physical value of the strange-quark mass the CFL phase has been found to be the ground-state 
at all potentially relevant densities. 

In this paper we will report on the results when incorporating neutrality conditions. As the CFL phase is neutral in 
the chiral limit due to the residual symmetry, we do not expect it to be influenced strongly. In contrast, the 2SC phase 
for two flavors is far from bein g ele ctrically neutral. Therefore, as known from investigations in Nanibu-Jona-Lasinio 
(NJL) type models fl§L fisl. [l6l. TlTj . the region of the 2SC phase is strongly affected by neutrality constraints. When 
varying the strange-quark mass the corresponding phenomenology can be, in a strongly simplified manner, sketched 
as follows: In the chiral limit, all quarks have the same Fermi momentum, and due to Luttinger's theorem, also have 
CO ' the same density. Increasing the strange-quark mass the Fermi momentum of the strange quarks is lowered according 
, to pf,s ~ M ~ -^s To first approximation the positive electrical charge of the up-quarks is then compensated by 
■ additional down-quarks with Fermi momentum pp.d ~ M + /2^. This mechanism is related to a partial symmetry 
QQ ' breaking of the residual symmetry in the CFL phase after the formation of Cooper pairs. Eventually, in some 
. channels ungapped q uasiparticle excitations appear, and accordingly these phases are referred to as the gapless CFL 
(gCFL) phase jll, [ll, Ho, HH. At some stage pairing of down- and strange-quarks is no longer favorable, and all 
remaining Cooper pairs involve an up-quark, which is characterizing the uSC phase. After this, the pairing of up- and 
strange-quarks might break, and one potentially may find in an intermediate regime the 2SC phase. Finally, all Fermi 
surfaces are far separated and the ground-state is no longer a spin-zero superfluid. Ungapped uSC and 2SC phases 
have not been observed in our calculations. It is worth mentioning that gapless phases are expected to be unstable, 
as some Meissner masses become imaginary [22, 23, 2^]. This finding led to a broad discussion of the appearance of 
inhomogeneous phases [25|, [2^, |2 7| . [2 8l . Due to the employed truncations this will not be addressed in this work. As 
suggested by our previous work [14| we find however those scenarios to be irrelevant for realistic quark masses with 
gaps as large as we obtain. 

This paper is organized as follows: In section [TT] we briefly present the theoretical framework which involves a 
fully self-consistent treatment of the quark propagator in the CFL-like phase including neutrality constraints. The 
role of chemical potentials especially in the CFL phase will be discussed in section IIIII In section IIVI we discuss the 
numerical results for the CFL-like phase at non-vanishing strange-quark masses. Finally, we summarize and conclude 
in section |Vl The color-flavor basis is given in an Appendix. 
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II. THE LANDAU-GAUGE QUARK PROPAGATOR AT NON- VANISHING CHEMICAL POTENTIAL 
A. The truncated quark Dyson-Schwinger equation at non-vanishing chemical potential 

In Refs. [3, [2§| a closed truncated Dyson-Schwinger equation (DSE) for the Landau-gauge quark propagator at 
non- vanishing chemical potential has been solved in the chiral limit and for finite strange-quark masses, respectively. 
In the present work we follow the same scheme and notations. In the Nambu-Gor'kov basis (see Ref. [J for a review) 
the normal quark DSE is coupled to an equation for the gap functions. To describe their structure we define 

^ ( ^°o''^ s^ip) = ^cs^{-prc ) ' 
^^p) - [r+ip) s-{p)^^cs+{~pfc ) ' 

with So being the 'bare' and S the full Nambu-Gor'kov propagator. The Nambu-Gor'kov self-energy is denoted by 
S. As the Euclidean action is real we furthermore have to impose $-(p) = 74$+(p)t74. The inverse 'bare' quark 
propagator in the presence of a static, isotropic and homogeneous 'external' gluon field with non-vanishing time 
component A4 is given by 

1 _ f -ip-j-i{p4 + ip.+ ^gA4)j4, + m \ 

Here m is the mass matrix and p the matrix implementing chemical potentials associated to global symmetries. The 
parameterization and determination of p and will be discussed below. The rcnormalization constants Z2 and Zip 
also appear in the DSE for the quark propagator 

s-'{p) = Z2So\p) + Zifj:{p) (5) 

which explicitly reads 



T± = -Zif[z2S^ VZifS^) (6) 
Employing the truncation of Ref. [2^ , the equations for the self-energy and the gap function then result in 



" ^p^ ^ -^"i (2^^-^^^ ('^^^''^'^ [-k^rm ^ i^^Tm ] ' 

where k = p — q. Here projectors transverse and longitudinal to the medium have been introduced. The functions G 
and F describe the corresponding medium modifications of the gluon propagator, see Eq. (21) of Ref. [2§|. They can 
be calculated once the coupling is given. Therefore the only input for the quark DSE are the running coupling as{k'^) 
and the renormalized current quark masses. 

As the objective of the presented work is to show that the CFL phase is the physical ground-state for realistic 
strange-quark masses, we use the most conservative choice for a.s(fc^) in our approach. This is the running coupling 
as{k'^) determined in DSE studies of the Yang-Mills sector [13, HH, [s^l . As is detailed in Ref. [13 it underestimates 
chiral symmetry breaking significantly in the Abelian approximation. In Ref. [2^ this coupling lead to the smallest 
critical strange-quark masses. However, we would like to emphasize that in contrast to NJL-tye models, our approach 
is much less sensitive to the choice of the coupling, an effect which can be traced back to the inclusion of the medium 
polarization. 

The rcnormalization constants are determined in the (chirally broken) vacuum. Due to the vertex construction 
employed, the quark-gluon vertex rcnormalization constant. Zip, cancels in the resulting renormalized equations. For 
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each flavor, we determine the quark wave-function renormaUzation constant, Z2, and the renormahzation constant 
relating the unrenormalized quark mass mo,q(A^) at an ultraviolet cutoff A to the renormalized mass mq(v) via 

mo,9(A2) = Z„,{u'^,K')mq{u), (10) 

by requiring 

^t^P)\p2=u^ = -ii> + mq{v) (11) 

at a renormahzation scale v. This corresponds to a momentum-subtraction [MOM) scheme, which results in some- 
what smaller numerical values for the quark current masses at a given renormahzation scale (usually taken to be 
2 GeV) . We simply ignore here the difference between MS and MOM masses because the effect is of the order of ten 
percent (when calculated within perturbation theory) and is thus much smaller than the uncertainty in the value of 
the current masses. 

It turns out that, as expected, the mass dependence of the quark wave function renormahzation constant, Z2, is 
negligible as long as the renormahzation scale is much larger than the mass. Therefore, we simply drop this dependence 
and Z2 is determined once and for all in the chiral limit. To keep the number of parameters as small as possible, we 
set the up- and down-quark masses to zero and vary the strange-quark current mass only. 

In the following we will restrict to spatially isotropic phases. In order to solve the DSE of the quark propagator it 
is advantageous to consider their color-flavor structure first. To get a self-consistent solution we choose suitable sets 
of matrices {Pi} and {Mi\ in color-flavor space, such that 

= ^T.^tij>)P^^ (12) 

Zjif ^ 

I 

^^^P) = -^J2'^tip)M^, (13) 

where we have introduced the renormalization-point independent component functions Tif{p) and (f)'^{p), which are 
matrix-valued in Dirac space. Full self-consistency is guaranteed in case a basis of all allowed matrices is considered. 
The dimensionality of this basis in a given phase depends on the residual symmetry in color-flavor space. For the 
CFL phase this will be detailed below, and the basis used is cxplicitely given in the Appendix. 
The Dirac structure of the self-energies in an even-parity phase can be parameterized by [ssl ] 

E+(p) = -z^I]+,(p)-*^pS+,(p) + S+,(p)-*74^Si,,;(p), (14) 
^tiP) = (74^0X.b)+74<^J,,(p)+05,,(p)+^0i,,(p))75, (15) 

where p = p/\p\, f = p • ^, ijjp = ujpj^ and ujp = ip4 + p,. Thus we finally need to solve a coupled set of integral 
equations for the energy- and momentum-dependent functions ^^bcd iip) ^^^'^ 't'ABCD iip)- This is done numerically. 

B. Color neutrality 

As has been discussed in [1^ [13, [H, [sH] , we need to allow for constant values of A4, whose DSE reduces to 

= 0. (16) 
As j44 is anti-Hcrmitian within our conventions, we define by 

j:t^aX'' -^^9A, (17) 
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the effective color chemical potentials. The color chemical potentials are then adjusted to obtain color neutrality. 
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C. Electrical neutrality and /^-equilibrium 

Electrical neutrality within QCD alone is meaningless. We need to include the electro-weak interaction and its 
particle content. The conserved electrical charge is then adjusted by its Lagrange multiplier /iq and the chemical 
potentials for the quarks, including the color chemical potentials, are given by 

flab,i] = yU SabSij + ^ ^J.d>^ib^ij + fJ-Q SabQij , (18) 
d 

where a,b = 1, 2, 3 denote color and i,j = 1, 2, 3 flavor indices. The charges of the quark flavors are encoded by the 
matrix Q = diagy (|, — ^, — ^) = ^r^ + -^^'t^- Thus the 'bare' inverse quark propagator takes the form 

S-^ip) = (-'P^^^-'^P^ + 'f'^^^ + ^ V ° -"T, ^ V (19) 

Compared to the strong interaction, we can consider leptons as non-interacting particles. Only the electrons will be 
relevant, as the relevant charge chemical potentials will not or only slightly exceed the mass of other charged leptons. 
With charge —1, their charge chemical potential is pe ~ ^l-i-Q- This also means 

= ^J■s = M« + ^J■e (20) 

for each color, which is usually referred to as /3-equilibrium. We consider the electrons as massless and their electrical 
charge density given by p^i = -^Hq. Electrical neutrality then enforces 

= 0. (21) 



D. Parameterization of the CFL phase 



A finite strange-quark mass leads to a partial symmetry breaking of the CFL symmetry in the chiral limit through 
the mass matrix m = \rns (l — -v/Sr^). If in addition /ig ^ due to the neutrality constraint, we can directly 
conclude from its generator Q — + '^/^'^^ that 

SUcM'i) SU,+v{2)®Uc+v{l) 

C/c+y(l)®C/e+y(l), (22) 

where the residual symmetry is generated by — A^" and rg — A^. Those form a Cartan subalgebra of the Lie 
algebra of SUc+v{i), i-e. a maximum set of commuting matrices. The most general choice of a basis {Pi\ and {Mi\, 
constructed as in Ref. [sj for a less complicated ansatz, is given in the Appendix. It consists of 15 matrices for {Pi} 
and {Mi}, respectively. 

To preserve the residual Uc+v{^) ® C^c+y(l) symmetry, we are only allowed to vary the color chemical potentials 
/U3 and /xg. The question is, whether this is enough to fulfill the requirement p"(a;) = for a = 1, ... ,8. To clarify 
this, we define the color-charge density matrix p by 



(23) 



This 3x3 matrix is symmetric and can be interpreted as the matrix of color-charges in the basis s}. With its 

help the condition in Eq. (I16p is reduced to 



Tr,(pA°) = 0. (24) 
As the Gell-Mann matrices form an orthogonal basis, we conclude that we need to adjust {/ia} to get 

P = fl, (25) 
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with p then being the quark number density. 

Varying fj,^ and fx^, we will surely achieve the requirement in Eq. (j24p if the matrix is diagonal. This is the case 
for all NJL-type investigations [l^, [3, [IB, [i3|- However, this is a further truncation on the self-consistency and 
enforced by hand. In general the matrix p is only symmetric, if we also consider the color-flavor structures P4, . . . , Pg. 
Therefore all pa for a = 1, ... ,8 need to be taken into account and the residual symmetry gets completely broken. 
The same arguments also hold if pg = 0, but p^ ^ or ps ^ 0. 

We conclude, that a finite strange quark mass induces in general backgrounds fields that break the CFL symmetry 
completely and the CFL phase can no longer be defined by a residual symmetry. It can only be defined by a continuity 
argument for the ground-state as a function of TOs that is CFL symmetric for = 0. 

Instead of introducing further color chemical potentials aside from p^ and /ig, we will estimate those as described 
below and choose p^ and /ig such that p'^ = and p^ = 0. This approach is similar to the truncation being used in 
NJL-type investigation and we could otherwise not constrain {Pi} and {Mi}. The ansatz becomes self-consistent in 
the 2SC and uSC phase. 



E. Estimating pa not in the Cartan subalgebra 



Varying p^ and ps only, we adjust the 3 x 3-dimensional real and positive charge-density matrix p to have equal 
diagonal elements and we define Tr (p) = -^Pp- The off-diagonal elements of p are strongly suppressed. After 
diagonalization we obtain p = diag^(/5i, p2, P3) = pD and, approximating the system by a free gas, we estimate 
Pi = "^Pi ■ For ^ fr^e gas we would therefore need to vary p^ and p^ in the new basis by A//3 = ^ {p2 — pi ) and 
Aps = ^~^{'^P3 ~ Pi ^ P2) in order to obtain a neutral phase. This will be used as an estimate of pa for a 7^ 3, 8 after 
the transformation into the old color basis. 

The eigenvalues of p can be estimated in an expansion in pij <§; pn for i ^ j and are given by pn and pn ± dp, 
with Sp = V P12P21 + P23P32 + P13P31- The corresponding Fermi momenta pi are approximately given hy pp and 
Pi^(l±(5/5/3/5ii). Therefore, after ordering the Fermi momenta, we have: Ap^ « ppSp/ipa and Apg, « 0. As for a free 
gas we also suspect the chemical potentials to transform like p = pD with DD^ = 1. Thus \\Ap\\ = \\p — pl\\ = 
\\p — pl\\ = ||A/i|| for any matrix norm and we obtain 

Pa < \\Ap\\ « 4^PF, a ^3,8. (26) 
opii 

We will come back to this estimate when discussing our numerical results for the chemical potentials. 



III. ELECTRONS IN THE CFL PHASE AND GAPLESS PAIRING 



A. Are there electrons in the CFL phase? 

It has been argued that the color- neutral CFL phase is automatically electrically neutral [38j. Therefore pq = 
and no electrons are allowed in the phase, which would have important consequences. The same result is also found 
in self-consistent NJL-type investigations [isl. [l7|. Since we find deviations, we first describe the reasoning in those 
models: 

Apart from the SUc{3) (E) S'[/\/(3), we can also consider the electrical Uq{1) generated by Q. In the dynamical 
symmetry breakdown SUc{3) (X) SUv{3) Uc+vi^) ^ J7c+v(l) also Uq{1) gets broken, however the phase is still 
symmetric under the so-called Uq{1) symmetry that is generated by 

Q = Q-lx^- J-A« 
2 2V3 

= i (r3 - A^-) + ^ (r« - n • (27) 

In the basis ((r, u), (g, d), (&, s), (r, d), (g, u), (r, s), (&, u), (g, s), (&, d)) used in the Appendix, the matrix Q is diagonal 
and the quasiparticles carry the charges (0,0,0,-1,1,-1,1,0,0), respectively. From this it has been concluded 
that Pq ^ in the fully gapped CFL phase, which then is a Q-insulator [l^]. The quark contribution to the 
thermodynamical potential Pq[T, p, p^, ps, pq] would then be invariant imder 

Pq[T,p,P3,P8,PQ] ^ Pq[T,p,ps-^p^,ps- ^^Pq,pq+ Pq]. (28) 
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Neglecting the leptons, we would therefore have degenerate ground-states under variation of ^g. Among those the 
real ground-state among those is then chosen by the minimum of the electronic (or leptonic) contribution to the 
thermodynamic potential, Pei[T', M3) Msj Mq] ^ TTK^f^Q- would therefore conclude jig = 0, which means that no 
electrons are allowed in the system. 

The whole argument is therefore based on the assumption pg — 0. Considering the form of the charge Q, we can 
concentrate on the separate pairing of {(r, d), {g, u)} and {(r, s), (6, u)}, respectively, as only those carry non- vanishing 
Q-charge (see Appendix). The requirement pg = is therefore equivalent to the statement, that two pairing fermion 
species with different chemical potentials have the same density in a fully gapped phase. We will now show that this is 
true for energy-independent gap functions, which are used in NJL-type models and emphasize that energy-dependent 
gap functions, as in our framework, do not require this. 

For the case of two fermion species a, b we neglect normal self-energies and consider 



n s- 



St J 



/ -if + /ia74 -75A* 

~ii> + tibli -75 A* 

75 A -ii>-lJ.ali 

\ 75A -ii> - pbli 



(29) 



giving 



[i-ajhl^ - I A 



2 + Mb/a74 



(P4 - i^^hjaY ^V^ ' 



(30) 



With the definition of Da/5 = (^4 + ^[i'alb)^ + and LOa/b = Pi + ij^a/bi the density Pa/b{x) of the fermions species 



and b turns out to be 



Pa/b{x) 



(^/I^Tr^ {slfbiPh*) 

d^p fdp, ^^Dl/a{^a/bDl/^+\A\^i^*b/a) 



For an energy independent gap function A, we can perform the energy integral and obtain for j/ia — ^b\ < 2|A| 



(31) 



Pa/bix) 



_fp_ 

(2^ 



'iP-fir + \A\ 



^{p+fif 



(32) 



where p, = ^ {fia + Hb)- The momentum integral is finite for a gap function that vanishing sufficiently fast and the 
result shows that for small differences in the chemical potential the densities of both particle species are the same. 
This has also been found in Ref. [s^l- For larger splittings, i.e., \p,a — lJ.b\ > 2|A|, gapless modes emerge (which will 
be discussed below). 

As we have a non-trivial energy dependence of the gap functions the above argument does no longer hold and we 
may find qualitative differences from NJL-type investigations. This may be directly verified by using an ansatz for 
the gap functions. The presence of electrons is therefore also connected to a finite width in the spectral function of 
the quasiparticles. The latter is included in our investigations and is discussed in Ref. (37j for the 2SC phase. 

We would like to clarify that a finite width of the quasiparticles near the Fermi surface may be an artefact of our 
current truncation, especially since we have not implemented the Meissner effect and the interaction via magnetic 
gluons is long-ranged. A closer investigation should be subject of future work. 



B. Gapless pairing 



The discussion of the last subsection sheds also light on gapless pairing. For illustration we will again use the simple 
parameterization in Eq. If the difference in the chemical potentials = (/i^ — Pb)/'^ exceeds the value |A|, the 

gap in the excitation spectrum vanishes and we find gapless modes even for A 7^ 0. Their dispersion relation are given 
by detS'+,^^ = 0, which are the zeros of the denominator of the integrand in Eq. (|3ip . For the CFL phase, gapless 
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Figure 1: Dependence of certain gap functions (^J ^ at the Fermi energy and for vanishing momentum (left) and of the (effective) 
cliemical potentials (right) on the renormalized strange-quark mass. Both dependencies are determined at jj, — 400 MeV. 



pairing might occur in some channels [39| and has been found as the preferred homogeneous ground-state in NJL-type 
investigations [isl. [l9l. [20l. [2l[ due to neutrality constraints. However, they are expected to be unstable as the Meissner 
masses become imaginarv |22l. [23l . US] ■ This finding led to a broad discussion of the appearance of inhomogeneous 
phases [2^, [23, US ■ We will postpone this subject to future investigations and consider here homogeneous phases 
only. The problem of imaginary Meissner masses in the gap equation has not been addressed yet and does not appear 
in our truncation either. 

For \5^\ > |A| gapless modes create a breached pairing region (40l |. For the occupation numbers implicitly given in 
Eq. (15^ and choosing /ia < /if,, we get 



n-a/bip) 



M + P T V (P + m) 



for fi - - |A|2 < p < /i -f - |A|2, 



2^(p + /2)^ + |AP 



(33) 



else. 



We see that in an interval of the Fermi momentum the occupation number is almost vanishing for one species and 
is almost unity for the other. The deviation of the occupation numbers from zero and unity in this interval for our 
simple parameterization is solely coming from the anti-quasiparticles and therefore of order A^/^x^ ^ \ xhe physical 
interpretation is that the quasiparticles in this 'breached pairing region' do not pair. 



IV. NUMERICAL RESULTS 



We will now present numerical results for the neutral CFL phase. The amount of data is overwhelming, as we 
calculate 4 x 15 dressing functions depending on energy and momentum for normal self-energy, gap function, normal 
propagator and anomalous propagator, respectively. However, we refer to Refs. [l3, for details and will only 
highlight new features coming from the neutrality constraints. 

In Fig. [1] we present the dependence of the gap functions on the renormalized strange-quark mass taken at the Fermi 
energy, i.e., for p4 = 0, and for vanishing momentum,. As the color-flavor structure of the pairing is quite involved, 
we refrain from a determination of the Fermi momenta as we have done in Ref. and focus on the gap functions 
at vanishing momentum. This is reasonable as the gap functions (/>J ^ for moderate chemical potentials are almost 
constant below the Fermi energy (see Refs. [l^ill^)- The gap functions shown belong to the tensor structures Pi, P2, 
P3, P4, Pe and P7 given in the Appendix and are labeled according to the quarks involved. As can be seen, the gap 
functions are weakly dependent on the renormalized strange-quark mass up to values of ^ 135 MeV. At this value, 
the pairing of green strange-quarks with blue down-quarks becomes ungapped as discussed in section 1111 Bl This can 
also be observed in the corresponding occupation numbers (see Fig. [2]). From the value 135 MeV onwards, the 
pairing gap of strange- and down-quarks 'melts' and above ~ 180 MeV we are in the uSC phase. This phase exists up 
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Figure 2: Occupation numbers of green strange-quarks (solid) and blue down-quarks (dashed) for various values of the renor- 
malized strange-quark mass at a chemical potential of jj. — 400 MeV (left) and the critical renormalized strange-quark masses 
as a function of the chemical potential together with the range stated by the particle data group [SJ] (shaded area) (right). 



to a value of ~ 190 MeV and in the window of ^ 190 — 200 McV we find the 2SC phase. Above this value no spin-zero 
pairing takes place anymore. 

For numerical reasons we have not been able to compute the pressure difference between the phases. However, 
the transitions of CFL to gCFL phase, and from there to the uSC phase, are likely to be at most of second order 
because the self-energy and gap functions change continuously. In addition, if the intrinsic symmetry is not altered, 
the transitions might even be crossovers. The transitions of uSC to 2SC phase, and of the 2SC to the unpaired phase 
look like first-order transitions. Here we also refer to our earlier observation, see Ref. [ij], that in these phases it is 
much more difficult to find an energetically disfavored solution. (As we solve the quark propagator DSE, being the 
variation of the CJT action, by iteration in the high-dimensional space of 'discretized' dressing functions, the domain 
of attraction for the global minimum is as usual strongly dominating.) More important, we want emphasize that 
those transitions are not in the physically relevant range of the strange-quark mass. 

On the rhs of Fig. [1] we present the electron chemical potential fie, the color chemical potentials ^3 and as well as 
an upper bound for the chemical potentials not belonging to the Cartan subalgebra being labeled M!{3,8}- It becomes 
obvious that electrical neutrality puts the strongest constraint on the phase structure and even exceeds the size of 
the largest gaps in the gCFL, uSC and 2SC phase. It is remarkable to note that /ig is non-vanishing in the gapped 
CFL phase as discussed in section IIIII This can also be seen in Fig. [21 which shows (in this case for green strange- 
and blue down-quarks) that the occupation numbers of pairing quarks in the gapped phase need not be identical as 
concluded for energy independent gap functions in Eq. (j33p . The bound for the chemical potentials not belonging to 
the Cartan subalgebra is of the order as the chemical /ia and /is . This shows that all color chemical potentials are 
much smaller than the electron chemical potential fig- Thus we estimate that the treatment of /is and /is alone is not 
much better justified than neglecting all color chemical potentials. 

In Fig. [2] the occupation numbers for the ungapped pairing of green strange- and blue down-quarks at a chemical 
potential of /i = 400 MeV and for various renormalized strange-quark masses are presented. Qualitatively those show 
a breached pairing region as found in Eq. ([33]) for a simple model study. Due to the interaction, these occupation 
numbers in the breached pairing region are neither almost vanishing nor close to unity, but more or less close to the 
occupation numbers in the unpaired phase (see also Ref. [29j). 

Last but not least, in Fig. [2|we present the critical values of the renormalized strange-quark mass, i.e., the values 
of the renormalized strange-quark mass which separate the different phases, as a function of the quark chemical 
potential. As compared to our analysis given in Ref. [l3| the transition from the CFL phase is slightly shifted towards 
higher values of the strange-quark mass when imposing neutrality constraints. This is indeed the anticipated result: 
The self-energy and gap functions in the CFL phase are only weakly modified by enforcing neutrality constraints, 
whereas these functions are strongly affected in the phases at large strange-quark masses. In agreement with the 
general arguments given in the introduction, towards larger strange-quark masses we first find the uSC phase followed 
by the 2SC phase. Above some critical value of the renormalized strange-quark mass none of the superfluid phases 
within the class of our ansatze is energetically favored anymore. This is the case because the Fermi surfaces of the 
different flavors are to far sep arated to allow for a spin-zero pairing. In this region one expects pairing of quarks 
with the same flavor only [U |4^ , see also Ref. [3] foi' a treatment in the Dyson-Schwinger approach. We also find 
a gCFL phase, which is supposed to be unstable (cf. the discussion given in the introduction). For the physical 
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strange-quark mass, given as a range of values by the particle data group [4l|, this becomes only relevant for the 
stated upper limit of the strange-quark mass and quark chemical potentials < 370MeV. For those values NJL-type 
calculations [13, E, [H, typically find the chirally broken phase. (In order to achieve a conservative estimate for 
color-flavor unlocking this does not yet take place in the presently used truncation due to our approximation for the 
medium polarization.) In addition, it is worth mentioning that solutions of the Bcthe-Salpeter equation within the 
DSE approach in the MOM scheme favor small values for the physical strange quark-current mass [45|, a fact which is 
also related to the difference in rcnormalization schemes. Taking into account all these facts wc conclude that among 
all possible color-superconducting phases only the CFL phase is realized for physically relevant strange-quark masses 
and chemical potentials. 



V. CONCLUSIONS 



As could have been anticipated from our previous investigation of color-flavor unlocking, the CFL phase stays the 
ground-state at vanishing temperatures and above the critical chemical potential for the chiral phase transition also 
when imposing neutrality constraints. The underlying reason for this result, which deviates from the one obtained 
from NJL-type models, is the inclusion of medium polarization effects. The latter play an important role in presence of 
a Fermi surface and lead to smaller dynamical mass generation due to damping and screening [l^ . As a consequence a 
stronger explicit symmetry breaking by the bare strange-quark mass as present in nature would be needed to separate 
the Fermi surfaces and to 'unlock' the CFL phase. To reinforce our conclusions we have used a coupling, which is 
rather weak and not able to give physical values for the chiral condensate and pion decay constant in the vacuum and 
which thus should favor color-flavor unlocking. 

Furthermore, it has been clarified that for a non- vanishing strange-quark mass the CFL phase cannot be defined by 
a residual global symmetry, and thus should better be named color-fiavor-locked-like phase. This can be traced back 
to the fact that non- vanishing chemical potentials /is, /is or induce static gluonic background fields which break 
the residual symmetry. This effect is, however, estimated to be small. 

A surprising feature is the appearance of electrons in the CFL phase which has been expected not to take place. 
Their presence is allowed due to the energy-dependence of self-energy and gap functions which give the quasiparticles 
a considerable width but may an artefact of the employed truncation. 
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Appendix: THE COLOR-FLAVOR STRUCTURE OF THE CFL PHASE 



For the neutral CFL phase with two degenerate quarks and non- vanishing strange-quark mass, wc generalize the 
ansatz of Rcfs. ([O, [13]) by choosing the matrices 



/ 5ii 5 a (5i6 

<^j5 5i2 5t7 
' Si9 Si3 



P,. = 



10 



'Jil3 



il4 



\ 



(A.l) 



il5 
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Si6 








Si8 


Si9 


Si3 



A/,; 



JilO 



Sill 



'Jil2 



s. 



13 
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(A.2) 



This is a complete basis for a residual Uc+v{^) ^ Uc+v{^) symmetry (see Eq. (HU). The basis is defined by 

{{r, u), {g, d), {b, s), {r, d), {g, u), (r, s), (6, u), {g, s), (6, d)}, 
with r, g, b denoting the color and u, d, s the flavor of the quarks. 
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